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ABSTRACT 


The  problem  considered  is  that  of  minimum  time  intercep¬ 
tion,  subject  to  fuel  constraints,  of  a  nonnaneuvering  target 
by  an  interceptor.  The  motion  of  the  interceptor  is  con¬ 
trolled  by  application  of  a  thrust  vector  (with  bounded  magni¬ 
tude)  and  by  the  control  of  the  thrusting  angle  relative  to, 

say,  the  interceptor  velocity  vector.  This  problem  is  solved 

2 

using  the  minimum  principle  of  Pontryagin  .  It  is  shown  that 
the  time-optimal  interception  policy  is  (a)  to  use  maximum 
thrust  in  the  beginning  of  the  control  interval  and  (b)  to 
use  a  constant  thrust  angle  during  the  thrusting  interval. 

The  optimal  thrust  angle  is  determined  analytically  as  a 
function  of  the  position  error,  the  velocities,  and  of  the  re¬ 
maining  fuel. 
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On  Minimum  Time  Exoatraospheric  Interception 
with  Fuel  Constraints :  A  Closed  Loop  Solution 

I  THE  EQUATIONS  OF  MOTION 

The  interception  problem  considered  takes  place  outside  the  atmosphere. 
Only  the  planar  version  of  the  problem  is  considered.  The  geometry  of  the 
problem  is  illustrated  in  Fig.  1. 

Consider  a  point  target  T.  Let  x^t)  and  y,p(t)  denote  the  horizontal 
and  vertical  coordinates,  respectively,  of  the  target  at  time  t.  Let  V^^Xt) 
and  V^j(t)  denote  the  target  horizontal  and  vertical  velocities.  Then,  the 
target  motion  is  defined  by  the  four  differential  equations 


ip(t)  =  V^Ct) 
yT(t)  =  v^Ct) 
v^t)  -  0 


(1.1) 


where  g  is  the  acceleration  of  gravity. 

Next  consider  the  interceptor  I.  Let  x^(t)  and  y^(t)  denote  its  position 
coordinates  in  the  x-y  plane  and  let  vxj(t)  and-  V^^t)  its  velocity  components 
and  let  m(t)  denote  its  mass.  Suppose  that  a  thrust  can  be  generated  by  the 
burning  of  fuel;  let  c  denote  the  exit  velocity  of  the  propellant  (assumed 
constant).  Let  P(t)  denote  the  rate-of-flow  of  fuel  of  the  propellant.  The 
thrust  f(t)  is  then  defined  by 

f(t)  =  c0(t)  (1.2) 


Let  ty(t)  denote  the  thrust  angle  with  respect  to  the  horizontal. 

Using  these  variables,  the  motion  of  the  interceptor  is  governed  by  the 
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following  set  of  differential  equations 
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Fig.  1.  The  geometry  of  the  interception  problem. 
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The  differential  equations  for  the  interception  problem  can  be  simplified  by 
defining  the  following  position  and  velocity  error  variables 

ex(t)  =  x^t)  -  ^,(t) 
e  (t)  =  y^t)  -  yT(t) 

A  M) 

vx(t)  £  VxI(t)  -  VxT(t) 

Vy(t)  £  VyI(t)  -  VyT(t) 

From  Eqs.  (l.l),  (1.3),  and  (1.4)  one  can  verify  that  the  error  variables  and 
the  mass  satisfy  the  differential  equations 

*x(t)  -  vx(t) 
ey(t)  =  vy(t) 

vx(t)  B  coe  t(t)  (1*5) 

vy(t)  -  ^  sin  *(t) 

itt) 

The  five -dimensional  vector 


3 


=  x(t)  (1.6) 


shall  serve  as  the  state  vector  for  the  optimization  problem. 

There  sire  two  constraints  that  will  be  imposed.  The  first  constraint  is 
that  the  thrust  that  is  generated  is  bounded  by  the  relation 

0  <  f(t)  <  F  for  311  t  (1-7) 

The  second  constraint  is  that  the  mass  of  the  available  propellant  is  limited. 

Thus,  if  mg  is  the  mass  of  the  interceptor  without  fuel  and  mQ  is  the  mass  of 

the  interceptor  fully  fueled,  then  the  mass  of  the  fuel  is  mQ  —  m  . 

The  problem  of  interception  in  minimum  time  can  now  be  posed  as  follows : 

Problem  1.1  Given  the  system  described  by  the  differential  equations 

(1.5) .  Suppose  that  at  the  initial  time  t  =  0  the  values  of  e  (0),  e  (0), 

X  y 

vx(0),  Vy(0),  and  m(0)  =  mQ  are  known.  Then,  determine 

(a)  the  thrust  f(t)  (0  <  f(t)  <  F),  and 

(b)  the  thrust  angle  \|/(t) 

so  that 


e  (T  )  =  e  (T  )  =  0 
x  y 

(1.8) 

m(T  )  >  m 

—  e 

(1.9) 

where  T  is  the  minimum  possible  interception  time. 

In  the  following  section,  the  minimum  principle'*'  will  be  applied  to  de¬ 
termine  the  necessary  conditions  for  optimality.  The  necessary  condition  will 
then  be  used  to  determine  the  optimal  thrust  and  the  optimal  thrust  angle. 


e  (t) 
xv 

e  (t) 

y 

V  (t) 
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V  (t) 

y 
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II  THE  NECESSARY  CONDITIONS  FOR  OPTIMALITY 


In  this  section,  several  relations  will  be  developed  which  characterize 
the  optimal  thrust,  thrust  angle  and  the  time-optimal  trajectory.  These  rela¬ 
tions  are  obtained  by  a  straightforward  application  of  the  minimum  principle 
to  the  Problem  1.1. 

Let  p^(t),  p^(t),  p^(t),  p^(t),  and  p,-(t)  denote  the  costate  variables 

(Lagrange  multipliers)  associated  with  the  state  variables  e^(t),  ey(t), 

v  (t),  v  (t),  and  m(t),  respectively.  Since  Problem  1.1  is  a  time-optimal 
x  y 

problem,  the  Hamiltonian  function  is  given  by 


H  =  1  +  ex(t)p1(t)  +  ey(t)p2(t)  +  vx(t)p3(t)  +  vy(t)p^(t)  +  m(-t  )p^(t )  (2. 1 ) 
which  upon  substitution  of  Eq.  (1.5)  yields 

H  =  1  +  vx(t)P1(t)  +  vy(t)p2(t)  +  P3(t)  cos  \|r(t) 

+  Sln  Pjfa)  (2.2) 

The  costate  variables  satisfy  the  differential  equations 


Px(t) 

P2(t) 

P3(t) 

i4(t) 

P5(t) 


Sh  n 

^7tT=  0 

Sh  n 

?yty  0 
"  sritT  -  - 

SH  /x  n 

■  STT tj  ■  -  ^ 

*y 

—  ‘SnT?T  =  ^  [p3(t)cos  f(t)  +  p^(t)sin  t(t)) 


(2.3a) 

(2.3b) 

(2.3c) 

(2.;d) 

(2.3e) 
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fran  which  one  immediately  deduces  that 


p^(t)  =  (constant) 
p2(t)  =  (constant) 

Ps(t)  = 

=  *4  -  «2t 


(2.4) 


The  posed  boundary  conditions  at  the  terminal  time  T  and  the  transversality 
conditions  yield 


ex(T*)  =  ey(T*)  =  0 

p3(t*)  =  p4(T*)  =  0 

(0  if  m(T*)  >  m 
P5(T  )  =  \  * 

?  (?  if  m(T  )  =  me  (or  f(T  )  =  0) 


(2.5) 

(2.6) 

(2.7) 


"ft 

Since  the  terminal  time  T  is  free  and  since  the  Hamiltonian  is  not  an  explicit 
function  of  time,  then  along  the  optimal  trajectory 


H  =  0  for  all  t  €  [0,T*]  (2.8) 

Let  \|f  (t)  and  f  (t)  denote  the  optimal  thrust  angle  and  thrust,  respectively; 
according  to  the  minimum  principle^-  they  must  absolutely  minimize  the 
Hamiltonian.  Noting  that 

sfftj  —  p3(t)  sln  +  p4(t)  cos  (2'9) 

"  1(1}  P3W  C0E  -  If!}  V*)  Bln  *<*>  (2.2.0) 
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then  it  is  necessary  that 


dH  I 


=  0 


(2.11) 


Mt  y 


>  0 


t(t)=4  (t) 


(2.12) 


From  Eqs.  (2.8)  and  (2.UL)  one  deduces  that 


tan  V*(t)  = 


p^tT 


(2.13) 


and  from  Eqs.  (2.10)  and  (2.12)  that 

P3(t)  cos  \|r*(t)  +  p^(t)  sin  V*(t)  <  0  (2.1  +  ) 

Equation  (2. 13)  yields  the  optimal  thrust-angle  in  terms  of  the  costate 
variables.  It  remains  to  determine  the  relationship  satisfied  by  the  thrust 
f  (t).  Towards  this  end,  define  the  switching  function  s(t)  by 

,  a  Pq^)  Pk(t)  PcS(t) 

s(t)  =  cos  *(t)  +  5(tT  sin  *(t)  "  (2'15) 

so  that  the  Hamiltonian  (2.2)  takes  the  form 

H  =  1  +  vx(t)  px(t)  +  vy(t)  p2(t)  +  f(t)  s(t)  (2.16) 


"K" 

The  requirement  that  f  (t) 

minimizes 

the  Hamiltonian  yields  the  relation 

F 

when 

s(t)  <  0 

f*(t)  = 

0 

when 

s(t )  >  0 

? 

when 

s(t)  =  0 
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Note  that  if  s(t)  =  0  over  a  finite  time -interval,  then  Eq.  (2.17c)  corres¬ 
ponds  to  the  singular  condition. 

This  set  of  equations  completes  the  list  of  relations  that  must  hold 
along  the  optimal  trajectory.  In  the  next  section,  these  relations  will  be 
used  to  deduce  additional  properties  of  the  optimal  thrust-angle  f  (t)  and  of 
the  optimal  thrust  f  (t). 


HI.  PROPERTIES  OF  i|r  (t)  AND  f  (t) 

-X* 

First,  it  will  be  shown  that  the  optimal  thrust  angle  i|r  (t)  must  be 
constant  throughout  the  time  of  thrusting.  This  conclusion  can  be  deduced 
by  the  following  argument. 

Recall  that  (see  Eq.  (2. A))  the  costate  variables  p^(t)  aad-  P4(t)  are 
given  by 


P3^  =  *3  “  Itit 
P4(t)  =  *4  -  *2* 


(3.1) 


Use  of  the  transversality  condition  (2.6),  i.e. 

p3(t*)  =  p4(t*)  =  0 


yields 


«2 


or 


Now  consider  the  relation  (2.13) 


(3.2) 


(3-3) 


(3*M 
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tan  \|f  (t)  = 


P^C*) 


(3-5) 


Ey  taking  time-derivatives  of  both  sides  of  Eq.  (3*5)  and  by  using  Eq.  (2.3) 
and  (3.1)  one  obtains 


It  tan  " 


~  P3(t)p^(t)  +  p3(t)p^(t) 


P3OO 


Pl ^ ~  ~  «2t)  ~  n2(^3  ~  n1t) 


p|(t) 


Pj(t) 


”ln4  **2*3  _  0 

P^t) 


(3.6) 


in  view  of  Eq.  (3.4).  This  implies  that 


tan  \|r  (t)  =  constant 


(3.7) 


and 


(t)  =  t  =  constant  (3.8) 

The  fact  that  \|r  is  constant  represents  useful  information.  It  cannot 
be  as  yet  evaluated  because  one  must  determine  the  thrust  f  (t).  For  this 
reason,  the  properties  of  f  (t)  will  be  determined.  It  will  be  shown  below 
that 

(a)  the  singular  condition  (2.17c)  does  not  occur, 

(b)  the  optimal  thrust  f*(t)  is  full-on  (f*(t)  »  F) 
at  the  beginning  of  the  time- interval  until  all 
the  fuel  is  consumed  at  which  time  it  is  turned 
off  (f*(t)  =  0). 
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The  fact  that  the  singular  condition  (2.17c)  cannot  occur  is  proved  by 
contradiction.  Let  ?  denote  a  subinterval  of  [0,T  ].  Suppose  that  the 
switching  function  s(t)  is  identically  zero  for  all  t  e  ?.  Thus,  suppose 
(see  Eq.  (2.I5)) 

P=(t) 


s(t)  = 


P3(t)  * 


m 


TO- 


cos  t  + 


p^oo 

m(t) 


sin  - ^ -  =  0  for  all  t  e  ? 

(3.9) 


where  we  have  used  the  fact  that  \|r  =  constant.  Since  the  Hamiltonian  must 
be  zero  along  the  optimal  trajectory  (see  Eq.  (2.8)),  substitution  of  (3*9) 
into  (2.16)  yields 


1  +  vx(t)  n-j.  +  v  (t)  n2  =  0  for  all  t  e  ? 


(3.10) 


Differentiation  of  (3.10 )  and  substitution  of  (1.5)  yields 


cos  V*  +  "2  sm  ■  0  for  all  t  e  ? 


Hence,  (3. 11)  implies 


(3.H) 


.  *  “1 
tan  \|r  =  —  — 

*2 


(3.12) 


But 


*  pjt)  n,  -  ic2t 

tan  \|r  =  -rrV  =  — - =r 

P3(t)  «3  -  «xt 


(3.13) 


Hence  (3.12)  and  (3.13)  imply  that 


[rtl  +  *2]  1  =  *1*3  +  *2*k 


(3.14) 


which  can  hold  if  and  only  if 
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nl  =  n2  =  0  (3.15) 

But  substitution  of  (3*15)  into  (3*10)  yields  0  «  1, obviously  a  contradiction. 
This  implies  that  the  hypothesis  s(t)  =  0  for  all  t  e  ?  is  false;  therefore, 
the  singular  condition  (2.17c)  cannot  occur.  This  in  turn  implies  that  the 
optimal  thrust  f  (t)  can  only  attain  its  full-on  value  f  (t)  =  F  or  its  cutoff 
value  f  (t)  =  0. 

The  next  topic  to  be  investigated  deals  with  the "shape" of  f*(t) .  Intu¬ 
itively,  one  would  expect  that  the  thrust  should  be  turned-on  as  soon  as 
possible  so  as  to  attain  the  maximum  velocity  possible.  This  intuitive 
feeling  can  indeed  be  verified  using  the  necessary  conditions. 

To  deduce  the  switching  nature  of  f  (t)  one  examines  the  switching  function 
s(t)  (see  Eq.  (3.9)).  Differentiation  of  s(t)  and  use  of  m(t)  =  -  and 

of  (2.3)  yields 


s(t)  = 


m2(t) 


f*(t) 


[p3(t)  cos  y*  +  Pj^t)  sin  i|r*] 


+  m(t)  [p3(t)  cos  V*  +  p^(t)  sin  i|r*] 


-  tPo(t)  cos  \|r*  +  p.  (t)  sin  >|r*] 

m2(t)c  3  4  (3.16) 


and  so,  Eq.  (3.16)  reduces  to 


s(t)  ■  [jtx  cos  \|r*  +  «2  sin  V*]  (3.17) 


We  compute  s(t)  to  deduce  that 


s(t) 


cos  v*  +  Itg  sin  Y*] 

m  (t)c 


(3.18) 
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From  Eqs.  (3.17 )  and  (3*18)  we  obtain  the  relation 


S(t) 


which  is  integrated  to  yield 


(3-19) 


s(t)  =  b(t)  exp 


/ 


4?  *> 

m(<J)c 


(3-20) 


Since  the  exponential  function  is  always  positive,  it  follows  that  s(t)  always 
has  the  same  sign  (positive  or  negative)  for  all  t  e  [0,T  ]. 

At  the  terminal  time  T  ,  the  switching  function  s(T  )  is  given  by 


s(T*)  =  — --  [pq(T*)  cos  \| *  +  pu(T*)  sin  i*] 
m(T  )  3  4 

Since  p^(T*)  =  p^(T*)  =  0  (see  Eq.  (2.6))  it  follows  that 


p5(t*) 


(3.21) 


* 

s(T  )  =  -  — 2 - 


(3.22) 


We  claim  that 


s(T*)  >  0 


(3.23) 


To  see  this  suppose  the  contrary,  i.e.,  that  s(T  )  <0.  Then,  Eq.  (3*22)  yields 
P^(T  )  >0.  But  let  us  recall  that 


P5(t)  =  [p3(-t)  cos  V*  +  fy(t)  sin  **]  50  (3# 


m  (t) 


2k) 
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in  view  of  the  necessary  condition  (2.14).  Hence,  if  p^(T*)  0,  then  Eq. 

(3*24)  implies  that 


P5(t)  >0  for  all  t  e  [0,T  ] 


(3.25) 


Since  the  switching  function  s(t)  is  given  hy 


s(t)  =  tP3(t)  cos  if*  +  p^  sin  **]  - 

then  the  necessary  condition  (2.14)  and  Eq.  (3.25)  yield 


Pc;(t) 


(3.26) 


s(t)  <»  0  for  all  t  e  [0,T  ] 


(3.27) 


But  this  means  (see  Eq.  (2.17a))  that  f  (t)  *  F  for  all  t  e  [0,T  ],  i.e.,  that 
full- thrust  is  applied  throughout  the  control  interval;  this  clearly  violates 
the  fuel  constraint  and  Eq.  (2.7).  Therefore,  s(T  )  <.0  cannot  occur  and,  so, 
Eq.  (3.23)  must  be  true. 

It  was  shown  above  that  s(T  )  >  0.  We  also  demonstrated  that  s(t)  always 
has  the  same  sign.  If  s(t)  <  0  for  all  t  e  [0,T  ],  then  s(t)  <  0  for  all 
t  e  [0,T  ],  which  means  that  f*(t)  =0  for  all  t  e  [0,T*];  but  if  no  thrust 
is  applied,  then  interception  in  general  is  not  possible.  It,  therefore 
follows  that 


s(t)  >  0  for  all  t  e  [0,T*] 


(3.28) 


and,  hence,  s(t)  is  a  monotonically  increasing  time -function.  This  means  that 
s(t)  <  0  during  the  beginning  of  the  time-interval  [0,T  ],  crosses  zero  at 
some  time  t  =  tQ  and  then  remains  positive.  This  establishes  the  fact  that 
one  must  thrust  in  the  beginning  until  all  the  fuel  is  exhausted  (at  tat) 

G 

and  then  coast,  i.e.. 
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(3.29) 


f*(t)  =  F  for  all  t  e  [0,tj 

f*(t)  =0  for  all  t  e  (t  ,T*] 


To  recapitulate:  It  has  been  proven  that  the  optimal  thrusting  program 

is  to  apply  full  thrust  as  soon  as  possible;  during  this  thrusting  period, 

the  angle  \|r(t)  of  the  thrust  must  remain  constant.  In  the  following  section 

this  angle  will  be  found  as  a  function  of  the  state  variables. 

The  thrust  cutoff- time  t  can  be  computed  because  at  t  =  t  all  the 

c  c 

available  fuel  has  been  exhausted.  Since  f(t)  =  F  for  t  e  [0,tc],  the 
differential  equation  for  the  mass  is 


m(t)  =  — 


Oil.-  l 

c  c 


(3.30) 


At  t  =  0,  m(t  ) 
'  o 


in  while  at  t  ss  t 
o 


,  m(t  )  =  m  , 
c  c7  e 


Thus, 


ni(tc) 


(3-31) 


and,  hence,  the  cutoff  time  t^  is  given  by 


t 

c 


[mQ  -  me]  c 
F 


The  coasting  period  is 


(3-32) 


IV.  A  CLOSED-FORM  EXIKESSION  FOR  THE  THRUSTING  ANGLE  V* 

Let  us  examine  the  situation  at  the  thrust  cutoff  time  t  .  Let  e  (t  ), 

c  x'  c  7 

e  (t  ),  v  (t  ),  v  (t  )  and  m(t  )  «  m  denote  the  values  of  the  five  state 
y  c  x  c  y  c  ce  ^ 

variables  at  In  order  to  accomplish  interception  at  t  «  T  ,  one  must  have 

ex(T*)  =  ey(T*)  =0  (4.1) 
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Since  there  is  no  thrusting  for  t  e  [tc,T  ],  it  follows  that  the  velocity- 


vector  at  t 


|^x|^c  jjmust  point  toward  the  origin  of  the  ex  —  plane  (see 


Fig.  2).  This  implies  that  at  t  =  t  the  following  relations  must  hold 

c 


e  (t  )  v  (t  )  a  e  (t  )  v  (t  ) 
x'  c'  yv  c'  y'  c'  x'  c' 


(1*.2) 


ex(tc>  vx<tc)  '  0 

W  W  '  0 


(*.3) 


These  relations  together  with  m(t  )  =  m  define  a  surface  S  in  the  five- 

0  6  C 

dimensional  state  space.  Clearly,  if  the  state  "belongs  to  S  ,  then  inter- 

•K*  ^ 

ception  is  possible.  In  fact,  the  time  T  —  t  can  be  evaluated  analytical  - 
and  is  given  by 


T  -  t  = 
c 


e  (t  ) 
xv  c7 

V  (t  ) 
XN  c ' 


ey(tc} 


(4.4) 


It  should  be  now  clear  how  one  can  determine  the  thrust  angle  \Jr  prior 

to  the  cutoff  time  t  .  Let  t  e  [0,t  ]  and  let  e  (t),  e  (t),  v  (t),  v  (t),  and 

c  c  x  y  x  y  ^ 

m(t)  be  the  values  of  the  five  state  variables  at  that  time.  Since  F  and  t 
are  constant  from  t  to  tQ,  the  differential  equations  (1.5 )  are 


*,<*>  -  \<*> 


k  <t)  = 
y 

V  (t)  = 

x'  ' 

V  (t)  = 

y 


V  (t) 

F 

mCtJ 

F 

m(t) 

F 


cos  \|r 

■> 

sin  ^ 


(4.5) 


m(t)  =-  i 


These  equations  can  be  integrated  from  t  to  t  to  yield: 
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Fig.  2.  Relation  of  velocities  and  errors  at  cutoff 
to  guarantee  interception. 
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m(tc)  =  me  -  m(t)  -  |  (tc  -  t) 


(4.6) 


v  (t  ) 
x'  c' 


"x(t)  +f  dT  =  Vx(t)  +  p  008 


dr 


t  »(t)  -  §  (T-t) 


vx(t)  +  F  cos  y*  .  |  log 


All 

F 


m(t)  -  £  (tc 


(4.7) 


Using  (4.6)  one  finds  that 


v  (t  )  *»  v  (t)  +  c  log  (slii 
Xs  c'  x'  '  V  me 


.  COS  t 


(4.8) 


In  an  identical  manner  one  obtains 


W  -  vy(t)  +  °  ^  l2^) 


sin  y 


(4.9) 


Finally 


t  t 

c(tc)  =  ex(t)  +  J  vx(t)  dr  >  ex(t)  +  J  jvx(t)  +  c  log  ^  --  j  cos  y*jdr 

"fc  c 


=  ex(t)  +  vx(t)[tc  -  t]  +  cos 


♦*..  i0**  r— ^ — •] 

t  U(t)  -  |(T-t)J 


dT 

(4.10) 


Evaluation  of  the  integral  and  use  of  (4.6)  yields 


e  (t  )  «  e  (t)  +  v  (t) [t  -  t]  +  c [t  -  t]  log  m(t)  cos  y*  +  c[t  -  t]  cos  y* 

A  C  A  A  v  V  V# 


+  me  log  me  •  cos  y*  -  m(t)  log  m(t)  •  cos  y* 


(4.11) 
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But  from  (4.6) 


tc  -  t  =  j;[m(t)  -  me] 


(4.12) 


Therefore,  e  (t  )  can  he  simplified  to 
X  c 


ex(tc}  =  ex(t)  +  f  [m(t)  "  me]  vx(t)  +  T  [me  log  +  m(t)  "  me]c°8  ^ 


2 

c 

F 


(^.13) 


In  a  similar  manner  one  obtains 


ey(tc)  »  ey(t)  +  |  [m(t)  -  me]  vy(t)  + 


T  [me  108  5(fy  +  m(t)  -  me]8in  * 


(4.l4) 

Substituting  Eqs.  (4.14),  (4.13),  (4.8),  and  (4.9)  Into  Eq.  (4.2)  one  obtains 
the  equation 

[ex(t)  +  |[m(t)  -  me]  vx(t)  +  y.  q(t)  cos  +  c  log  [^^]*  sin  t*j 


[< 


=  I  ey(t)  +  |[m(t)  -  me]  vy(t)  +  y.  q(t)  sin  t*J*[vx(t)  +  c  log|"^-|«  cos 


e 

(**■.15) 


where  q(t)  is  defined  for  the  sake  of  notational  simplicity  by 


q(  t )  =  -  me  log^2^-)-*-  - 


m 


(4.16) 


By  performing  the  indicated  multiplications,  Eq.  (4.15)  reduces  to 


B(t)ex(t)  +  r(t)vx(t) 


sin  f*  -  ^6(t)ey(t)  +  r(t)vy(t) 


cos  if 


=  ey(t)vx(t)  -  ex(t)vy(t) 


(4.17) 
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where  8(t)  and  ir(t)  are  defined  by 
8(t)  S  c  log  (  Sill) 


/.  X  A  c: 
r(t)  =  — 


m(t)  log  -  m(t) 


+  m 


(it. 18) 


(!*.19) 


In  this  manner,  the  problem  reduces  to  the  solution  of  Eq.  (4.17)  for 


^  .  .fe  note  that  Eq.  (4.17)  is  of  the  form 


7T  Tv 

i^(t)  sin  ijr  +  a^(t)  cos  )Jr  =  b(t) 


(4.20) 


where 


a-L(t)  =  6(t)ex(t)  +  r(t)vx(t) 


a2(t)  =  -  |  8(t)ey(t)  +  r(t)vy(t)( 


b(t)  =  e  (t)v  (t)  -  e  (t)v  (t) 

y  x  Ay 


(4.21) 


(4.22 


(J+.23) 


But 


*  *  f  p  ~ 

a1(t)sin  +  a„,(t)cos  t  =  /a^t)  +  ag(t) 


sin 


*  -1  a2(t) 

*  +  tan  S^tT 


(2 


Therefore,  from  Eqs.  (4.24)  and  (4.20)  we  deduce  the  desired  solution 


=  sin 


-1 


b(t) 


yf 


-  tan 


■1  2 


ao(i) 


raf(t)  +  a*(t) 


qrtT 


(*1.25) 
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If  the  relation 


b(t) 


\l' 


a^(t)  +  a^t) 


<  1 


(4.26) 


holds,  then  interception  is  possible.  Otherwise,  the  interceptor  will  miss 
the  target. 

The  complete  structure  of  the  control  system  that  generates  the  optimal 

*  , 
thrust  angle  is  illustrated  in  block  diagram  form  in  Figs.  3  through  6  . 
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Fig.  3*  Generation  of  error  variables  [cf .  Eq.  (1.4)] 


Fig.  4.  Generation  of  auxiliary  time  variables  q(t),  6(t), 
and  r(t)  [cf.  Eqs .  (4.l6),  (4.l8),  (4.19)].  nig  is  the  empty 
interceptor  mass,  c  is  the  equivalent  exit  velocity,  and  F 
is  the  maximum  thrust. 
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Fig.  5 .  Generation  of  auxiliary  time  functions  an(t),  a2(t), 
and  b(t)  [cf.  Eqs .  (4.21),  (4.22),  (4.23)]  from  tne  error 
variables  (see  Fig.  3)  and  the  auxiliary  variables  y(t)  and 
6(t)  (see  Fig.  4) . 
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RUN  OUT  OF  FUEL;  NO  NEED  TO  FIND  V' 


YES  =>  COMPUTE  THRUST  ANGLE  'f' 

SET  THRUST  TO  FULL  ON  f*(t)  =  F 


Fig.  6.  Generation  of  optimal  thrust  and  thrust  direction  angle  <Jj*. 
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